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We present a convenient analytical parametrization of the deuteron wave function calculated 
within dispersion approach as a discrete superposition of Yukawa-type functions, in both configura- 
tion and momentum spaces. 
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Recently in the paper [l[ it was shown that the 
deuteron tensor polarization component T2q(Q 2 ) pro- 
vides a crucial test of deuteron wave functions in the 
range of momentum transfers available in to-day experi- 
ments. The calculation [l[ shows that the most popular 
model wave functions do not give adequate description 
of T2o(Q 2 ) and are to be discriminate in favor of those 
obtained in the dispersion potentialless inverse scattering 
approach with no adjustable parameters Q, (see Q, too) 
and giving the best description. Some time ago this func- 
tion (MT - wave function) was used in the calculation of 
the neutron charge form factor Q ■ The results of calcula- 
tion (12 new points) are compatible with existing values 
of this form factor of other authors. A fit is obtained for 
the whole set (36 points) taking into account the data 
for the slope of the form factor at Q 2 = 0. These re- 
sults will be used in the neutrino scattering experiments 
in Fermilab ||. 

The aim of the present paper is to present a con- 
ventional algebraic parametrization of the deuteron MT- 
wave function calculated within dispersion approach as a 
discrete superposition of Yukawa-type terms. 

Let us remind briefly the main characteristic features 
of these wave functions obtained in the frame of the po- 
tentialless approach to the inverse scattering problem. 

The important feature of these wave functions is the 
fact that they are " almost model independent" : no form 
of NN interaction Hamiltonian is used. The MT wave 
functions are given by the dispersion type integral di- 
rectly in terms of the experimental scattering phases and 
the mixing parameter for NN scattering in the 3 Si — 3 Di 
channel. Regge-analysis of experimental data on NN 
scattering was used to describe the phase shifts at large 
energy. 

It is worth to notice that the MT wave functions were 
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obtained using quite general assumptions about analyti- 
cal properties of quantum amplitudes such as the validity 
of the Mandelstam representation for the deuteron elec- 
trodisintegration amplitude. These wave functions have 
no fitting parameters and can be altered only with the 
amelioration of the NN scattering phase analysis. 

Let us notice that the process of constructing of these 
wave functions is closely related to the equations ob- 
tained in the framework of the dispersion approach based 
on the analytic properties of the scattering amplitudes 
(| 0] (see also [|[ and especially the detailed version 
9(). In fact, this approach is a kind of dispersion tech- 
nique using integrals over composite-system masses. 

Let us emphasize that by construction the dispersion 
wave functions [2] differ principially from deuteron wave 
functions used in the conventional nuclear model (see, 
e.g., the review [l(|). So, for applications it is convenient 
to have MT-wave functions [2j in analytical form. 

Therefore, we present here a simple perametrization 
of the deuteron function as a superposition of Yukawa- 
type terms (that was introduced in Ref. [Ill for Paris 
potential; see also the fit in Ref. [l2} for CD-Bonn wave 
function). 

So, we consider the deuteron wave functions (pi(r) in 
the states with orbital momentum I = ,<po(r) = u(r) 
and I — 2 , ^{r) = w(r) The ansatz for the analytic 
versions of the r-space wave functions, denoted by u a (r) 
and w a {r), is 
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rrij = a + mo (j — 1) 



where the coefficients Cj , Dj , the maximal value of the 
index j and mo are defined by the condition of the best 
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FIG. 1: S-wave deuteron wave function. Solid line - MT @,0|, 
dushed -[H, dot-dushed - [El. 




FIG. 2: D-wave deuteron wave function. Legend is the same 
as in FigfT] 



fitting, a = ^/M Ed, M is nucleon mass, Ed is the binding 
energy of deuteron. 

These wave functions are normalized according to 



dr 



(u(r)f + (w(r)) 2 



1 



(2) 



The conventional boundary conditions at zero: 

u(r) ~ r , w(r) ~ r 3 , (3) 

lead to one condition for Cj and three constraints for Dj , 
namely: 



$> = o, E fl ^E^ = E3=°- ( 4 ) 
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Using the form (Q} it is easy to describe the standard 
behaviour of the deuteron wave functions at r — > oo. The 
asymptotics of S'-state is 



u(r) ~ A s e ° 
and the asymptotics of D state is 



rjAs 1 + — + —— 



(5) 



(G) 



TABLE I: Coefficients for the parametrized deuteron wave 
function calculated within a dispersion approach. The last Cj 
and last three Dj are to be computed from eq. (|4]|(n u =16, 
n w =12). 



0.87872995+00 
-0.50381047+00 

0.28787196+02 
-0.82301294+03 

0.12062383+05 
-0.10574260+06 

0.59534957+06 
-0.22627706+07 

0.59953379+07 
-0.11282284+08 

0.15181681+08 
-0.14519973+08 

0.96491938+07 
-0.42403857+07 

0.11092702+07 
eq. © 



D,(fm- 1/2 ) 

0.22245143-01 
-0.41548258+00 
-0.18618515+01 

0.21987598+02 
-0.16885413+03 

0.76001430+03 
-0.22287203+04 

0.43330023+04 
-0.54072021+04 
eq. © 
eq. © 
eq. © 



Here A$ and Ad = r\ As are the asymptotic S'-state and 
-D-state normalizations and r\ is the asymptotic " D / S 
state ratio" . In our calculation of MT-wave-function we 
use a=0.231625 far 1 . 

The Fourier-transforms of wave functions tpi (k) , I = 
0, 2 in the momentum representation in r-space are: 



k 2 dk 3l {kr)^{k) 



(7) 



where ji(kr) is the spherical Bessel function. 

The normalization condition for these wave functions 
is given by 



k 2 dk 



1 



The fits for momentum space wave functions, following 
from ((T|) and ([7]) are 



C, 



7 v J ; 



r 2 {k) 



E 



D, 



(fc 2 + m 2 ) 



(9) 



The calculated coefficients in the fits ([T]) , ^ are listed 
in the table 1. 

The asymptotics at r — > oo gives for the obtained fits 
of MT-wave functions the asymptotic "D/S state ratio" 



n 



0.02531511 



(10) 



The accuracy of this parametrization is illustrated by 
the magnitudes of the integrals: 



dr [u(r) - u a (r)Y 



1/2 



= 4.1-10" 



(11) 



3 



{roo \ 1/2 

J dr[w{r)-w a {r)] 2 \ = 2.2 • HT 3 . (12) 

So, we present a convenient analytical parametrization 
of the deuteron wave function calculated within disper- 



sion approach as a discrete superposition of Yukawa-type 
functions. This function is plotted on the Fig.l and Fig. 2. 
The wave functions [Tlj and fl2| are given for compari- 
son. 
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